Wigner function for squeezed negative binomial state and evolution of density operator for amplitude decay
1. IntroductionSqueezed states and their practical applications have been widely explored during the past several decades (see, e.g., the review articles[1,2] and references therein). This is mainly because the noises of squeezed states may be below the vacuum or ground-state noise level in a certain observable, which can be useful for interferometric detection and gravitational waves by means of improving the interferometric and spectroscopic measurements. Especially, two-mode squeezed states served as a new entangled resource is crucial in realizing some vital quantum information processing, such as quantum teleportation,[3,4] quantum computation,[5] and cryptography,[6] and can be readily generated by using experimentally controllable and available sets including phase shifters, beam splitters, and squeezers.[7] For another side, any system in the nature is not absolutely isolated, quantum noise usually occurs and causes the decoherence of system when this system is in the interaction with its surrounding environment. Hence, the environment-induced decoherence aspects in open squeezed optical fields have received great attention from researchers. For instance, the effect of amplitude decay, phase damping, and thermal noise on the decoherence of squeezed states has been analytically investigated in the pioneering works.[8–10]
In probability theory and statistics, the negative binomial distribution reads
, where
and s is a non-negative integer. Corresponding to this definition, there exists a pure state with a negative binomial form of its photon-number probability distribution, which is called the negative binomial state (NBS),[11] that is
where
is the Fock state obtained by repeatedly acting the photon-creation operator
on the vacuum state
in the Fock space. Using the identity
, we can rewrite
as
Indeed, the NBS is an intermediate state between a coherent state (
and
) and a quasi-thermal state (
s = 0 ). The reason for designating the state
for
s = 0 as a quasi-thermal state is that its photon number distribution is the same as that of a mixed thermal state, which is easily proved to be true using the compact form
, regarded as a multi-photon-subtracted thermal state via subtracting
s photons from a single-mode thermal field
with the mean number of photons
.
[12] This means, through subtracting
s photons, the photon number of
shows a negative binomial distribution. For the NBS and its derived states, various representations, properties, and quasiprobability distributions have been studied in some detail and their generation schemes in a number of nonlinear processes have been demonstrated.
[13–15] Especially lately, the applications of the NBS in quantum optics and its time-evolution in some typical decoherence channels including the amplitude decay channel (ADC),
[16] the diffusion process,
[17] and the laser channel
[18] have been investigated.
As a practical derived state, a squeezed negative binomial state (SNBS) was first introduced by Joshi and Lawande, and its second-order correlation function, Wigner function, Q functions, and photon number distribution have been investigated, and one result that needs to be noted is that the SNBS can be helpful for the transfer of quantum phase-space distributions from squeezed coherent states to squeezed number states.[19] In this paper, we shall introduce a new formula of SNBS, which is defined as an NBS
having undergone a two-mode squeezing transform, hence its density operator reads
where
is the single-mode squeezing operator with the squeezing
r.
[20–22] To go a step further, we focus on the theoretical investigation of Wigner function (WF) for SNBS and the evolution behavior of density operator in ADC by fully using the operator-ordering method
[23–26] and the thermal entangled state representation (TESR).
[27–30] Our motivations for focusing on this topic are explicated as follows. First, the operator-ordering method can help us prove that the WF for SNBS is expressed as the finite-sum representation determined by the non-negative integer
s, rather than the infinite-sum form as shown in Ref. [
19], which can more accurately describe the nonclassicality of SNBS in phase space (see Section
3 for details). Second, according to our best knowledge, the investigation of the evolution of SNBS in ADC by deriving the analytical results has not been previously discussed because of the complicated calculations resulting from the squeezing operation in SNBS, however TESR can help one to greatly simplify it and derive an accurate and concise result related to two operator Hermite polynomial functions.
2. Operator-ordering products of
To analyze the nonclassicality of SNBS from the aspects of phase space and decoherence evolution, we first use the Bogolubov transformation caused by the squeezed operator
to change SNBS into
where
, a weighted quantum superposition of squeezed number states
. Next, we derive the normal and antinormal ordering products of
. For this purpose, using the completeness relation of the coherent state
to derive the normal ordering of the operator
, that is
where the marked symbol
refers to normal ordering.
[31] Further, using the integral formula (17) in Ref. [
32] to perform the integration (
5) over
α within normal ordering, and operating this integral result on the vacuum state
and using the generating function of
n-order Hermite polynomial
, we have the compact form of the squeezed number state
, that is
which indicates that a squeezed number state can be viewed as a Hermite polynomial excitation on the squeezed vacuum state. Noting the normal ordering product of the vacuum projector,
i.e.,
, we directly obtain the normal ordering of the density operator
,
i.e.,
Using the summation formula concerning two single-variable Hermite polynomials in Eq. (10) in Ref. [
33], thus the normal ordering expansion of
is obtained as the compact exponential form
where we have set
,
, and
. Using a method similar to that used to obtain the result (20) in Ref. [
33], we can also prove that the normal ordering operator
is a bivariate distribution operator in phase space. Further, noting the operator identity
, equation (
8) can be rewritten as
which shows that the operator
is indeed a squeezed thermal state by comparing with the result (11) in Ref. [
33].
Further, substituting the normal ordering of
in Eq. (8) into the antinormal ordering formula[34] and using the mathematical integral formula,[35] we can obtain the antinormal ordering of
as the following form
where
denotes antinormal ordering
[36] and we have set
,
, and
.
3. Wigner function for SNBSThe Wigner distribution function is a useful tool which is used to investigate the nonclassical optical fields from the view point of phase space. In this section, we shall investigate the nonclassicality of the SNBS by analytically and numerically studying its phase-space WF distributions. For the given density operator
, its WF definition is given by
, where
is the coherent state representation of single-mode Wigner operator, i.e.,
.[37] Considering that
are coherent states,[38,39] hence we first introduce the normal operator identity of the operator
,[10] that is
and use the antinormal ordering of
in Eq. (
10), thus the state
reads
where
and
. Thus, inserting Eq. (
12) into the WF definition and using the completeness relation of the coherent state
and the generating function of
,
i.e.,
we therefore have
Further, using the integral formula[40,41] to complete the two integrations over β and z in the whole phase space, the WF
can be calculated as
where we have set
with
Further, using the exponential expansion of cross product term
, it then follows that
Noting the generation function of
in Eq. (
13) and its
l-order differential relation of
, we finally obtain
which is the analytical WF of the SNBS related to the square of the module of
-th single-variable Hermite polynomial. Obviously, when
s = 0,
becomes the WF for squeezed quasi-thermal field,
i.e.,
Mathematically, it is well known that two-variable Hermite polynomials
can be introduced either by its generating function
or by its power-series expansion
where
represents the transfer of the amplitude of number state in the time evolution process of a forced Harmonic oscillator,
[42] and can be useful for studying light propagation modes in quadratic graded-index media.
[43] Hence, a comparison of the right sides of Eqs. (
15) and (
18) directly leads to the result that
If we replace
,
, and
in Eq. (22), we therefore have
where
is a new two-variable special function whose generating function is
, a new and interesting result. Thus, the differential relation of the special function
can be obtained as
which is the same form as the well-known differential relation of
.
[29] Further, using the special function
, we find that the right side of Eq. (
22) can be represented as
, thus the analytical WF for SNBS can be simplified as
| |
To clearly observe the nonclassicality of SNBS, we plot its WF in (Re
) phase space for a fixed squeezing r and different values of the parameters s and γ as shown in Fig. 1. Obviously, the SNBS always shows squeezing effect in one of the quadratures for
as a feature of nonclassicality of this state. For
(e.g., s = 100) and
(or s = 0), the WF becomes the standard Gaussian distribution of squeezed coherent state (or squeezed quasi-thermal state) [see Figs. 1(a) and 1(d)], which is completely accordant with the analytical result outlined in Ref. [19]. Comparing Figs. 1(a), 1(b), and 1(c), we find that, with the increase of s, the WF gradually shows an inverted peak in the center of phase space, however it slowly loses the inverted peak as γ increases and finally becomes Gaussian when γ tends to one. Moreover, the larger value of s can lead to the smaller peak value of WF, but the variation of the peak value with γ is just a fully contrary trend. This result can easily be understood via analyzing the result in Eq. (2) since more and smaller number states are annihilated and the probabilities of finding number states in this field becomes small with increasing s (especially at center part in the phase space), however the larger γ causes the higher weights of the existing number states. Besides, for the limited case (
), the WF distribution seems basically unchanged with the increase of s, which is because the value of parameter s has a negligible effect on the weights of the existing number states.
4. Evolution of SNBS in ADCFor any open system, noise always arises in fundamental dynamics processes. For instance, the purely dissipative noise in ADC leads to the nonclassicality deterioration of the system. In this section, we shall analytically obtain the evolution formula of SNBS in ADC and investigate how does the amplitude dissipation affect the nonclassicality of SNBS.
4.1. Kraus operator solutions to master equation describing ADCIn the interaction picture, the density-operator evolution of the system in ADC can be characterized by[44–47]
where
κ is the decay rate of the system. Indeed, this process completely differs from the phase-decay process with a purely dephasing noise. To obtain the solutions of the above operator equation, we review the two-mode TESR describing quantum entanglement between system and heat reservoir, which is defined in enlarged Fock space as
where
is a displacement operator with the complex variable
η, the operators
and
are respectively the real and fictitious creation operators, which yield the commutative relations,
i.e.,
and
. Under the basis of the entangled state
, there exist the following extremely useful interchange relations between the real and fictitious modes,
i.e.,
,
, and
, which play a vital role of solving the master equation (
27) since it can help convert the equation for the density operator
in Eq. (
27) into the evolution equation for the state vector
, that is
whose formal solution is
Fruther, through calculating the scalar product
, we find that it is proportional to the inner product between an exponentially decaying entangled state
and the initial state vector
,
i.e.,
, and its proportionality coefficient is just a Gaussian decaying factor
,
. This point clearly reveals the decoherence behavior of the channel, namely, the amplitude decay. To take off the density operator
from
, we need to use the completeness relation
, thus after strict derivations (for the detailed formula derivations see Ref. [
27]), its infinitive Kraus operator-sum representation of
as the solutions of the master equation (
27) can be calculated as
where
is the Kraus operator of
. Hence, for a known initial state
, the evolved density operator
in ADC can be freely obtained from Eq. (
31). Also, using the operator identity
, we can prove the normalization condition
.
4.2. Density-operator evolutionNoting the marshalling sequence of the operators
,
, and
in Eq. (31), we need to introduce the anti-normal operator identity of the operator
,[32] that is
and use the normal ordering of
in Eq. (
8). Hence, the state
reads
Substituting Eqs. (
8) and (
33) into Eq. (
31) leads to the form
where we have set
and
. Inserting the completeness relations of coherent states
and
into Eq. (
34),
[38,39] and using the operator identity
and the generating function of
in Eq. (
13),
[48] we obtain
Using the integration technique within an operator-ordering product
[49–52] to perform the two integrals over
α and
β, then
becomes
where
can be represented as
since the normal ordering symbol within another normal ordering symbol can be canceled,
i.e.,
, and the parameters
h,
h1,
h2, and
h3 respectively are
Comparing the right sides of Eqs. (
15) and (
36), and using the new special function in Eq. (
24), we immediately obtain the analytical evolution of SNBS in ADC, that is
which means that the pure SNBS evolves into a new mixed state in ADC, related to two operator Hermite polynomial functions whose variables are the linear combinations of creation and annihilation operators within normal ordering. Especially, when
r = 0, thus
g = 1,
, and
, equation (
38) becomes the evolution of NBS in ADC,
i.e.,
with
, which shows that, after undergoing ADC, the resultant state is still an NBS, but the parameter
γ decays to
and
. It is in agreement with the result in Ref. [
16]. For the case of
s = 0,
becomes
which is the evolution of squeezed quasi-thermal field in ADC. In addition, when the decay time
, thus
T = 0,
h = 1,
,
, and
, so equation (
38) returns to SNBS in Eq. (
3). For the limited decay time,
i.e.,
,
, thus SNBS finally decays to vacuum.
5. ConclusionsIn summary, we have used TESR and the operator-ordering method to derive the analytical evolution law of SNBS in ADC. It is found that the analytical WF is the finite-dimensional sum over the square of the module of single-variable Hermite polynomials. Numerically, the WF for SNBS can exhibit a similar Gaussian distribution in both coherent and quasi-thermal limited cases, and the larger s can result in an inverted peak with the smaller peak value in the WF phase space but the larger γ plays an opposite role. For another case, after passing through ADC, the pure SNBS evolves into a new mixed state, whose normal ordering product of the density operator is related to two operator Hermite polynomials of the linear combinations of creation and annihilation operators, however the SNBS finally decays to vacuum as a result of amplitude decay when the decay time
is long enough. From the above results, we have clearly found that, in the study of the environment-induced decoherence aspects, the TESR and the operator-ordering theory provide a privileged tool to investigate the analytical evolution of density operator and characterize destruction level of nonclassicality. Besides, through deriving the WF for SNBS, we have found a new two-variable special function
and its generating function, with which the concise formulas of WF and density-operator evolution can be obtained.